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We calculate the ultrasonic absorption and the thermal conductivity in the superconducting state 
of UPts as functions of temperature and direction of propagation and polarization. Two leading 
candidates for the superconducting order parameter are considered: the E\g and E2u representations. 
Both can fit the data except for the ultrasonic absorption in the A phase. To do that, it is necessary 
to suppose that the system has only a single domain, and that must be chosen as the most favorable 
one. However, the E2u theory requires fine-tuning of parameters to fit the low temperature thermal 
conductivity. Thus, transport data favor the E\g theory. Measurements of the thermal conductivity 
as a function of pressure at low temperature could help to further distinguish the two theories. 



PACS numbers: 74.70.Tx, 74.20.Rp, 74.25. Ld 



I. INTRODUCTION 

UPt3 continues to pose a fundamental challenge to 
solid-state physics. It becomes superconducting at about 
Tc+ = 0.5if , entering the A phase, and has a second tran- 
sition at about Tc_ = QAAK entering the B phase. Its 
unusual phase diagram and thermodynamics show that 
it is unconventional. The precise nature of its supercon- 
ducting order parameter remains elusive, however. 

From the purely theoretical point of view, order param- 
eters (OPs) that transform according to a single represen- 
tation of the Z?g/j group and that take into account spin- 
orbit coupling are the most desirable. This includes the 
Eig [D and i?2u Hi theories. Because of the perception 
that these theories may not give sufficient phenomeno- 
logical flexibility, other pictures have been constructed. 
One of these is the mixed-representation {e.g., A-B) ||] 
picture. This suffers from the problem of fine tuning, 
since the coupling constants of the two representations 
must be nearly identical to account for the nearness of 
Tc+ and Tc- ■ Another picture is the extended-symmetry- 
group spin triplet theory Q. However, this requires zero 
spin-orbit coupling in an /-electron system, an artificial 
assumption. 

From the point of view of comparison to experiment, 
the most stringent constraints on theory come from de- 
tailed measurements of the thermal conductivity k |^ 
and the ultrasonic attenuation a |^ . Because of the di- 
rectional information inherent in these probes, they offer 
the best chance to distinguish between order parame- 
ters that may have only relatively subtle differences [Q. 
This expectation is reinforced by the fact that the ultra- 
sonic attenuation has been measured in the A phase, the 
most anisotropic of the zero-field phases. Previous cal- 
culations have indicated that the mixed representation 
and extended symmetry group scenarios are probably in- 
consistent with these measurements Accordingly we 
focus in this paper on the Eig and E2u theories. Our 
aim is to compare theory and experiment for k and a in 
a quantitative fashion. In so doing, we can also inves- 
tigate whether fine-tuning is required in either of these 



pictures. 

In the next section, we give the details of the calcula- 
tion. The following section compares the fits for the two 
models. The final section discusses how the models can 
be distinguished on the basis of pressure experiments and 
gives a brief summary. 



II. CALCULATION METHOD 
A. Bands 

The Fermi surface of UPts is highly complex. There 
are five or six sheets, some having a very complicated 
shape. However, the effective masses of the sheets have 
been very accurately determined for most directions by 
de Haas- van Alphen experiments These measure- 
ments show that the sheet contains about 80% of the 
density of states at the Fermi energy ep- This sheet is 
centered on the T point. We approximate it as an ellip- 
soid with a dispersion relation: 




(1) 



with rriz = 43.2me and rrix = rriy — n.-img. The 
small "eggs" of the Fa band at the zone boundary are 
ignored. We have investigated the effect of including 
other bands and find that the contributions are small. 
The anisotropic mass in the Fa band is very important, 
however. The results are quite different for an isotropic 
Fermi surface. 



B. Order Parameters 

The central part of the F3 band does not cross the 
Brillouin zone boundary. This means that nodes coming 
from the fact that gap functions should be invariant un- 
der translation by a reciprocal lattice vector do not need 
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to be considered. In this case, to good accuracy, one can 
use polynomial functions to describe the gap functions. 

Standard calculations for the Eig model use a singlet 
gap function 

A(k) = Ao(T) [Kk, + 1 5{T)Kky] , (2) 

while for the £'211 model conventional calculations use the 
triplet gap function 

d(k) = Ao(T)z [k,{kl ~ kl) + 2i6{T)k,k,,ky] , (3) 

where fc^ is the /^-component of k vector on the Fermi 
surface and z is the direction of the vector d(k) in pseu- 
dospin space. Both models involve the idea that the A 
phase involves only a single OP component: in the A 
phase, 5{T) = 0. Accordingly, S{T) is a temperature- 
dependent function that begins to grow continuously at 
the A-B phase boundary. 

A crucial point in the theoretical treatment, how- 
ever, is that the expressions in Eqs. (||) and (||) are 
not the complete OPs for these representations. Yip 
and Garg have shown that the complete gap func- 
tions for these phases are linear combinations of two 
pairs of functions for Eig rather than the single pair 
{kzkx , kzky) shown in Eq. (||). For i?2u, the gap function 
is a linear combination of six pairs rather than the single 
pair (zfcz(fc^ — ky) , zkzkxky). Writing the complete gap 
functions in the form appropriate for use in both the A 
and B phases, we have the formulas 

A(k) = Ao(T) ^ (k) + z^(T)/(")(k)] (4) 

n=l 

for Eig and 

d(k) = Ao(T) J2 «^"'[fi"Hk) + ^<5(T)f(")(k)] (5) 

for £'211. The functions /(k)and f (k) are defined in Ta- 
bles I and H, adapted from Ref . . In the Eig repre- 
sentation, the gap function is a linear combination of the 
two functions in the first column of Table |[ In the B 
phase, the functions in the second column also come in, 
with a relative phase of i, as shown in Eq. (|^). In the £211 
representation, the gap function is a linear combination 
of the six functions in the first column of Table ||. In the 
B phase, the functions in the second column also come in, 
with a relative phase of i, as shown in Eq. (0). Note that 



the resulting gap function can be multiplied by another 
completely symmetric (Aig) function without changing 
its symmetry. Nearly all calculations for Eig have set 
g^^-' = 0, while nearly all calculations for E2U have set 
— _ y^(4) _ y(5) _ y(6) „ f^YoY an exception 

see Ref. 

It is very important to understand what this special 
choice of g*-"^ or u^") implies. There are two levels of 



TABLE I: Basis functions for the Eig representation. 
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doing phenomenology for the OP. One is choosing a rep- 
resentation and working out the consequences. In this 
case the coefficients g*-"^ and u*^"^ must be treated as pa- 
rameters. They are continuous functions of all the many 
parameters that are in the microscopic Hamiltonian, as 
well as external variables such as pressure. Thus choos- 
ing a representation means choosing a set of Hamilto- 
nians of finite measure, i.e., to selecting an alternative 
that has a finite a priori probability. The other (more 
commonly used) level of phenomenology is to choose a 
specific basis function from Table | or Table |^ and work 
out the consequences. This means setting all but one of 
the 5^"^ or u^") equal to zero. It also means choosing 
a set of Hamiltonians of zero measure, which has zero 
a priori probability. In other words, choosing a specific 
basis function amounts to fine tuning. None of the g'"' 
or u^"^ can vanish by symmetry, since all symmetries of 
the system have already been taken into account in the 
group-theoretical decomposition. 

This does not necessarily invalidate previous work that 
has been done with a specific choice of basis function. If 
the comparison to experimental data is not very sensitive 
to variations in f/^"-* or then the fine tuning may 

be abandoned without harming agreement of theory and 
observation. A key goal of this paper is to investigate 
this sensitivity for the Eig and £211 theories. 

C. Transport Coefficients 

Ultrasonic attenuation experiments in UPta have been 
done in the hydrodynamic regime where the sound fre- 
quency Lo satisfies w < 1/ts , with Ts the electronic relax- 
ation time. The appropriate formula for sound propaga- 
tion along the direction q with polarization e ± g is: 



(T) oc / du 





dF{Lo) 




kf(e.k)^ 



(6) 
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TABLE II: Basis functions for the E2u representation. 



fr'(k) fr'(k) 

Ti — 1 kx X ky y ky X ~\~ kx y 

Til — 2 ( ky^kz z 2 kxkykz z 

Th — 3 {^kx ^kxk^^^ X ^^y^x^ y i,^x ^kxk.^^ y ~\~ ^^y^x^ 

Ti — 4 {^kx ^kxk,^^ X ~\~ i^ky ^kykj^^ y i^kx ^kxky^ y (^^y ^^y^x^ 

Tb — 5 ^kx ^k^ky ~\~ ky^kz z i^^kxky ^kxky^kz z 

Ti 6 i_k,j2 ~ \{)k^ky ~\~ ^kx ky 'jx ~\~ (^ky \Qkykj; ~\~ ^kykj^^y {kx ~ \Qkj^ky ~\~ ^kxky^y {^y lOA^^/c-p ~|~ ^kykj-^x 



Here F is the Fermi function and the angle brackets indi- 
cate a Fermi surface average. | Akp — > |dkP in the triplet 
case. The thermal conductivity in the ith direction when 



The prefactors are fitted to the values of a^i and kh at 

Our approach does not calculate the density of states 
self-consistently. This means that it is not very accurate 
at low temperatures, T < 0.2Tc+, but in this paper we 
are chiefly concerned with the behavior closer to Tc+. 

The temperature dependence of Ao(r) is not accu- 
rately known. However, one may deduce the power 
law at low temperature from the pattern of nodes, 
while near Tc we can use Ginzburg-Landau theory for 
the same purpose. For computational convenience, we 
take S{T) = [1 - (r/r,_)6]i/2 and Ao(T) = Ao(r = 

0) [1 - (r/rc+)="]^/^ . a; = 2 and X = 3 for point and hue 
nodes, respectively. We have verified that these functions 
are close to the Ginzburg-Landau result near Tc+ and 
they have the appropriate behavior at low temperature 
as well. In addition, this leads to a fit for the temperature 
dependence with only the single parameter Ao(T = 0). 
Finally, we follow Graf et al. Q in taking the empirical 
form for the relaxation time: ~ OmTTTc+{l+TyT^_^). 



D. Domains 

In the A phase there is an additional complication. 
6{T) = so that, as stated above, only the /x(k) (or 
£r(k)) functions are used. However, the theory requires 
that the superconducting OP is oriented by its coupling 
to the staggered antiferromagnetic moment Mg. There 
are three domains allowed by symmetry for M5, and neu- 
tron scattering experiments show that they are equally 
populated. The pattern of nodes rotates by tt/6 around 
the c-axis every time the wave passes through a domain 



a temperature gradient is imposed along the ith direction 
is computed using: 



(7) 

I 

wall. The gap function becomes, for example, /a;(7?.k) 
where T^k is the rotated vector. We shall take a coordi- 
nate system which is fixed in the laboratory: the propa- 
gation vector or polarization vectors refer to this frame. 

The observed a must result from an average over do- 
mains. In such a situation, response functions are usually 
calculated by effective medium theory, which may take 
various forms, depending on the approximations made 
[|l2|. However, the present case is simplified by the fact 
that a is the relatively small imaginary part of a large 
complex quantity, the sound propagation vector. For 
such a perturbative quantity, all averaging methods lead 
to the result that the observed global a should be the 
mean of the three local values. 



III. RESULTS 
A. Best fits 

We begin the comparison of theory and experiment 
in this subsection by adjusting all parameters in Eqs. (^) 
and (HI) to their optimum values for both representations, 
varying the mass anisotropy mz/rrix and choosing the 
most favorable domain of the A phase. Fairly impres- 
sive fits of data to experiment for both theories may be 
obtained in this way. 

The best fits for Eig and £'211 models for the ultrasound 
and the thermal conductivity are plotted against the data 
in Fig. |l|. The curves for k provide reasonable fits over 
the whole range of temperature. In particular, the overall 
anisotropy is well fit. Some discrepancies may be noted 
at sufficiently low temperature. The fits for a are also 
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FIG. 1: Best fits of Eig and i?2u models against tfie data 
[ p^ for both transverse sound attenuation (top frame) and 
tfiermal conductivity (bottom frame) . The fitting parameters 
are given in the text. 



FIG. 2: The figure shows the sensitivity of the fits to the 
cliange of the anisotropic mass ratio. In the bottom frame, the 
two curves (6 and c) are indistinguishable for the mz/nix — 1 
case. 



reasonable except that some minor discrepancies may be 
noted near Tc in the A phase. The best-fit parameters 
for Eig are: A(0) = 3Tc+, x = 2, and g^^) = 1.0 and 
g(^^ = 0.2, while the best-fit parameters for E2u are: 
A(0) — 4:.5Tc+, X — 2, and all the coefficients m^*^ vanish 
except u'-^-' — 1.0. 

There exists a major difference between the best fits 
of Eig and E2u models when anisotropic masses are con- 
sidered. For Eig, the anisotropic mass ratio is fixed as 
iTiz/'m-x = iTiz/my = 2.5. This agrees with the measured 
value 1^ for the dominant band. For E2u, in contrast, 
the best fits are acquired at rriz/mx = mz/niy = 0.33, in 
very poor agreement with the measured value. 

B. Anisotropic mass 

In order to give some idea of the sensitivity of the fits 
to the anisotropic mass ratio, we calculate a and k for 
E2u model using the same best-fit parameters except for 
the anisotropic mass ratio mz/mx- As seen in Fig. as 
the ratio is increased from 0.33 to its physical value of 2.5, 
the polarization anisotropy in a decreases sharply, while 



the anisotropy in k actually reverses. This indicates that 
the fit for the E2u model is somewhat questionable. 

Our fitting procedure differs from that of Graf et al. 
1^, |l3|. These authors used revised order parameters 

A(k) = Ao(r) [kzkx + i S{T) kzky] (1 + azfcf + a^kt) 

(8) 

for Eig and 

d(k) = Ao{T)z[kz{kl - kl) + 2i6{T)kzkxky]{l + a2kl + mkt) 

(9) 

for i?2ti, where the coefficients a2 and a4 are adjustable 
fitting parameters. Since the additional factor is com- 
pletely symmetric, this docs not change the representa- 
tion. As opposed to using the standard basis functions 
(a2 = 04 = 0), it is evident that large 02 and 04 will cor- 
respond to more electronic excitations around the line of 
nodes at k^ = 0. This in turn will enhance the in-plane 
thermal conductivity and sound attenuation relative to 
the out-of-plane ones. Roughly speaking, the larger 02 or 
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FIG. 3: The best fits of ultrasound attenuation data taking 
into account three equal-populated domain average. The best 
fits for one single domain (see also Fig. |^) are included for 
comparison. The curves for xz and {xz) are indistinguishable. 



a4 case for the revised model is physically analogous to 
having rtix , ruy > rriz (and hence larger in- plane density 
of states) in the standard model. This change in the gap 
function amounts to multiplying the relative contribu- 
tion of the in-plane excitations by a factor of 12 relative 
to what one would expect from the measured ratio of 
effective masses. 



C. Domain-averaged fits 

In order to fit the data in the A phase, it is neces- 
sary to choose the most favorable domain. This was also 
found by Graf et al. As mentioned above, this is not 
consistent with the assumptions of the theory. The idea 
that the order parameter is oriented by M5 is a central 
assumption for both pictures, since it accounts for the 
fact that Hc2 is isotropic in the basal plane. Thus do- 
main averaging is not optional: it is necessary for the 
consistency of the theory. 

If the absorption is properly averaged over all three 
domains, then we obtain the best fit for the two theories 
shown in Fig. ||. In the A phase, the fits are substantially 
worse than those one gets by choosing the best domain, 
as was done for Fig. |l] and in Ref. [||. 

The reason for the difference in the quality of fit in 
the A phase is mainly due to the node patterns. In the 
B phase, for both Eig and E2U models, there are lines 
of nodes along the equator and point nodes at the poles 
of the Fermi surface. In this phase, the nodal pattern 
is symmetric about the c-axis and thus the real nodal 
structures are virtually the same in the three different do- 
mains. In the A phase, in contrast, the two point nodes at 
the pole have connected to form line nodes in the fc^; = 



plane for Eig and the = iky planes for i?2u model (in 
the domain defined by the basis functions of Tables | and 
. In this case, the nodal structure is asymmetric about 
the c-axis. Depending on the orientation of the node pat- 
tern relative to the direction of polarization/propagation, 
the three domains contribute differently. As a matter of 
fact, as compared to an isotropic state, transverse sound 
with q II i and e || y (denoted as xy) will cause fewer 
excitations in the 7r/6 or 7r/3 rotated domains. It is this 
that causes the (xy) (the angle brackets indicate a three 
equal-populated domains average) attenuation to drop 
significantly in the A phase relative to the xy attenuation. 
In contrast, the best fits for the xz and (xz) attenuations 
are almost indistinguishable. These considerations apply 
equally to the Eig and E2U models. 

D. Fine tuning 

An important criterion for the acceptance of a theory is 
its robustness. If small changes in the parameters would 
spoil the agreement with the data, the theory is not ac- 
ceptable since it lacks explanatory power. Accordingly, 
we investigate the sensitivity of the OP's transport char- 
acteristics to changes in its shape. We shall focus on the 
B phase for reasons that will become clear below. 

For Eig, varying the shape is relatively straightfor- 
ward. There are only two parameters g^^^ and g^'^\ In 
the B phase with g^^^ = there are lines of nodes along 
the equator and point nodes at the poles. Changing the 
ratio of g'^^^ and g^^-' does not affect this nodal pattern. A 
larger value of g^^^ merely gives additional modulation to 
the order parameter in the azimuthal direction. We show 
some representative curves for Kc/ki, for different values 
of the ratio in Fig. ^(a). There is very little change in 
the goodness of fit as the ratio is varied. Defining a chi- 
square parameter for the fits as a function of g^^^ / g^^\ we 
see from its approximate constancy that the Eig theory 
is robust for this particular quantity. 

The E2U order parameter behaves very differently. For 
starters, the shape space u^^\ is six-dimensional. 

This precludes an extensive search. However, mapping 
of a two-dimensional subspace will suffice to make the 
point about parameter sensitivity. The main difference 
between E2U and Eig is clear from the nodal pattern. 
Taking into account u'^^'> in addition to the usual u'^^\ 
when u^'^-' = there is a line of nodes at the equator and 
point nodes at the poles in the B phase. When u'^^^ is 
finite, however, the line of nodes disappears. This is a 
consequence of Blount's theorem |Q, which states that 
triplet states in a spin-orbit coupled system do not have 
lines of nodes. In Fig. ^(b), we show some representative 
curves for Kc/ Kb for different values of the ratio m^^^ /u^^^ . 
The curves are very sensitive to the ratio, in particular 
at low temperatures. The chi-square parameter varies 
rapidly as a function of u^^^/u^^^. We conclude that the 
E2u theory is not robust. Its virtues in describing the 
data arise from a fine-tuning of the parameters. 
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FIG. 4: (a) The ratio of Kc/tt for different ratios of g'^^ /g'^^ 
for the -Big modeL (b) The ratio of /tc/wb for different ratios 
of u^^^ /vF^^ for the E2u model. 

IV. CONCLUSION 

Order parameters belonging to the two-dimensional 
representations Eig and E2u remain as the main candi- 
dates for the superconducting state of UPts. A difficulty 
that remains for both these theories is the explanation 
of the ultrasonic absorption in the A phase. The the- 
ories require that the observed absorption be the result 
of an average over three domains in which the order pa- 
rameter has three orientations. This averaging is delete- 
rious to the agreement between theory and experiment. 



This may indicate that there is something wrong with 
the idea that the magnetization orients the superconduc- 
tivity, or it may indicate that there is an additional scat- 
tering mechanism, such as scattering from domain walls, 
which partially cancels the drop in absorption due to gap 
formation. 

The theory has an additional severe problem at 
low temperatures. A line of nodes at the equator of 
the Fermi surface is rcqiiircd in order to give agreement 
with experiment for the thermal conductivity. In the six- 
dimensional parameter space of this theory, this occurs 
only in a two-dimensional subspacc. The theory therefore 
requires fine-tuning to explain the data; all the parame- 
ters in the microscopic interaction must conspire to give 
the desired result. It has been suggested that strong spin- 
orbit coupling can lock the d- vector along the z-direction. 
However, d does not have the interpretation of zero spin 
projection, as it does in ^He. Taking into account the 
spin-orbit coupling reduces the order parameter space to 
six dimensions for triplet E2, not one or two. 

Since the shape parameters m^"^ depend on all the de- 
tails of the microscopic Hamiltonian, the low temperature 
power laws that depend on having a line of nodes arc un- 
stable. By contrast, the line of nodes at the equator in the 
E\g theory is forced by symmetry. Application of pres- 
sure, for example would lift these nodes in an E2U gap- 
The pressure scale that governs superconducting phe- 
nomena may be estimated as Tc+{dTc+ / dP)~^ ~ Ikbar. 
One would expect very strong pressure dependence of the 
thermal conductivity on this scale at low temperatures 
for E2u but not for Eig. This could serve to further dis- 
tinguish the two theories. 

The data presently give preference to the Eig theory 
since it does not require fine tuning to explain the data. 
The discrepancy of theory and experiment for ultrasonic 
absorption in the A phase has not been resolved, and 
may require new ideas. 
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